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AN ANALYSIS OF LONGITUDINAL STABILITY IN POWER-OFF FLIGHT 
WITH CHARTS FOR USE IN DESIGN 

By Chablbs H. Zimmerman 



SUMMARY 

This report presents a discussion of longitvdinal 
stability in gliding flight together with a series of charts 
with which the stability characteristics of any airplane 
may be readily estimated. 

The first portion of this report is intended, primarily, 
for students of the subject. The relationships governing 
stability characteristics are derived from equations of 
equilibrium referred to moving axes that are tangent 
and perpendicular to the instantaneous flight path. It 
is shown that instability of the_ motion can arise only 
through an increase of linear and angular momentum in 
the system during one complete cycle. The interaction 
of events leading to increase or decrease of momentum 
during a cycle is explained in detail. The construction 
of charts showing the effects of the nondimensional 
dO dO 

parameters 0 L , 0 D , -j^> —m t , and —fim„ upon 

the stability characteristics is explained and the effects 
of the more important of the aerodynamic and mass 
characteristics of the airplane, as revealed by the charts, 
are discussed. 

The latter portion of the paper is devoted to a series of 
Jfi related charts with which the dynamic stability of any 
airplane in power-off flight may be readily estimated. 
The use of the charts is explained in detail so that 
reference to the earlier discussion is unnecessary. 

INTRODUCTION 

The longitudinal stability of aircraft has received 
very extensive and exhaustive treatment by able 
writers (see references and bibliography), but the 
classical treatment of the subject has been rather 
difficult for those not familiar with higher mathematics. 
The study reported herein was undertaken with the 
purpose of making more understandable the mathe- 
matical treatment and of preparing a method of esti- 
mating stability characteristics that would be suffi- 
ciently accurate and rapid to appeal to practical 
designers. 

The section preceding the group of charts for deter- 
mining stability characteristics in power-off flight is 
devoted to a derivation, in relatively simple terms, 
of the mathematical relationships and to a discussion 
of the formulas. The portion following the group of 
charts consists of an explanation of the method of 
using them. It is not necessary to read the first por- 
tion in order to use the charts with satisfactory results. 



All symbols not given in the report cover are denned 
where used and are also listed in the appendix. 

L ANALYSIS AND DISCUSSION 
DERIVATION OP MATHEMATICAL FORMULAS 

Definition of stability characteristics. — The sta- 
bility characteristics of an airplane are those qualities 
which define the nature of the motion after a deviation 
from an initial condition of equilibrium. The motion 
may be periodic, consisting of a series of oscillations 
having a certain period and rate of increase or decrease 
in amplitude, or aperiodic with a certain rate of return 
toward or deviation from the' eqmlibrium position. 
In many stable airplanes the return to a condition of 
nonoscillating equilibrium is spoken of as aperiodic or 
"dead beat" when it is essentially oscillatory in 
character but very heavily damped. 

Fundamental concepts and assumptions. — The forces 
and moments deter minin g the motion of the airplane 
are of two kinds: (1 ) Aerodynamic forces and moments 
created by movement of the lifting and control sur- 
faces relative to the surrounding air; (2) mass forces 
and moments arising from the weight and acceleration, 
angular as well as linear, of the airplane. The funda- 
mental basis of the discussion presented in this report 
is that at all times there exists a state of equilibrium 
between the mass forces and moments and the aero- 
dynamic forces and moments. 

A complete treatment of the stability of airplanes 
would be extremely lengthy and very complex. Cer- 
tain assumptions have therefore been made. As the 
motion of an airplane is three dimensional, it is to be 
expected that any treatment of the subject will be 
incomplete if it neglects certain of the components of 
the motion. Fortunately, conventional airplanes are 
symmetrical (within limits here applicable) with 
respect to the plane that includes the fuselage axis and 
is perpendicular to the span axis. It is obvious that a 
longitudinal motion having no component of linear 
velocity perpendicular to that plane or no component 
of angular velocity about any axis lying in that plane 
cannot introduce asymmetric forces or moments. 
Such motion can therefore be treated as an independent 
phenomenon. 

The longitudinal-stabihty characteristics will neces- 
sarily be affected by any deflection of the lifting or 
control surfaces. The influence of wing elasticity and 
of free longitudinal control will not be considered in 
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the primary analysis because such a consideration 
would complicate the relationships and obscure the 
fundamental principles. 

It is assumed that forces and moments acting upon 
the wing and the horizontal tail surfaces vary as the 
square of the air speed and the first power of the angle 
of attack of the individual surfaces and that they are 
not affected by the rate of change of either the air 
speed or the angle of attack. Also the forces upon 
the lifting surfaces are assumed not to be affected by 
the rate of rotation of those surfaces (reference 1). 

Each of the foregoing assumptions necessarily in- 
volves a certain degree of approximation but they are 
confirmed by comparison between measured and cal- 
culated values of stability characteristics (reference 2 
and unpublished data) and are justified by the simpli- 
fication of the relationships they permit. 

Equations of equilibrium. — As has been previously 
stated, the course of the airplane in flight is deterrnined 
by the conditions necessary to maintain equilibrium 
between mass and aerodynamic forces and moments 
at all times. In steady flight the equfiibrium may be 
expressed by the equations (see fig. 1): 




Fiqtoe 1.— Angular and vectorial relationships In flight, power off. 

W sin y+^pVSO^O a 
W cos y-^pV^SO^O b 
± P V*ScC m =0 



(1) 



where (la) refers to forces tangent to the instantane- 
ous flight path, (lb) to forces perpendicular to the 
instantaneous flight path in the plane of symmetry, 
and (lc) to moments about an axis through the center 
of gravity and perpendicular to the plane of symmetry. 
After displacement from the steady-flight condition 
the equations of equflibrium read: 

W sin (y+Ay)+±p(V+AV)*S(C B +AO D )=-m^ , 
PFcos(y+A 7 )— ^(V+AV) 1 S(C l + AC i )= -mV& b 



±p(V+AWSc(C n +AC m )=mk y $$ c 



(2) 



where 



dV 



AO„ may be written as 



acceleration tangent to flight path. 

ydy centrifugal acceleration normal to the flight 
dt' path. 

cPB angular acceleration of airplane about the 
dt 1 ' lateral axis. 

Since the effects of angular velocity and of acceleration 
upon the forces are neglected, AC D may be written as 

Aa^and A&asAa^; 

dC dO 

A a ^ g " + A g~gj^ where Ag=g, the angular velocity in 

pitch, since g is zero in the original condition; and 
AV, Ay, Aa, and g are small quantities by assumption. 
Terms involving products of two or more small quan- 
tities will be neglected. 
Also 

sin (7+A7)=sin y cos A7+COS y sin A7 
=sin 7+A7 cos 7 

and 

cos (7+A7)=cos 7 cos A7— sin 7 sin A7 
=cos 7 — A7 sin 7 

Then: 



W sin y+WAy cos y+±pV s SC D + P V&C D AV+±pV t g^Aa=-m^ 



W cos y-WAy sin y~^pV s SGr.- P VSG l AV— \ P V i ^Aa=-mvijt 



Tt 



^ P V 2 ScO„+pVScO n AV+^pV 2 Sc^Aa+^pV 1 Sc^i=mk Y ^ 



(3) 



Subtracting (1) from (3) 



WA7 COS y+pVSC D AV+±pV 2 g^Aa=-m^ 



'dt 



-WAy sin y-pVSCr,AV-l P V*tf^A*= -^9i 



dt 



P VScO m AV+^Sc^+^S^Aa= m k Y ^ 



(4) 
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From equations (1) 



Wcos y=± P V*SC L 
Wsmy=~pV s SC D 

c m =o 

Therefore, 

l P V i SC D A 7 -pVSC I AV-±pV>S^Aa=-mV d £ b 
\ P VS^ + \pV>S^Aa=rnk Y ^ c 



(5) 



Dividing (5a) and (5b) by ± pVS and (5c) by pSV 



AV 



0 D Ay-2G L - 



V 

AV dC L 



dO DA 0 m dt 

-d^ La - 2 pSv V 



.a« i m <kl 

V da WV dt 



^ b 



Let 



m 
7SV~~ 



aV =av> 



=r (after Glauert, reference 3) 



(6) 



T 
dv 

dt 
V" 



dV 
' dt 



—CiAy—ZCiAV'- 



dO DA 9 dV 
da at 



GoAy^OzAV ~^Aa= -2& 
da dt 



2 Vc -dJ^ + 2 Vc ^a- Aa==Tk ^ 



JO, 



(7) 



Equations (7a) and (7b) are nondimensional and the 
only variables are A7, AV, and Aa. In order to re- 
duce (7c) to a comparable form it is necessary to 

rewrite which is the only constant in the 
eg 

equations not in readily usable form, and to express 
0 in terms of y and a. 

The greater portion of the change in pitching mo- 
ment produced by rotation occurs because the rota- 
tion changes the angle of attack of the tail surfaces. 
A positive rotation g causes the tail to move downward 
with a velocity, relative to the center of gravity, equal 
to gl where I is the distance from the center of gravity 
to the mean quarter-chord point of the horizontal 
surfaces. The tangent of the change in angle of 
attack at the tail is tan Aa,=gl/V which may be re- 
placed by Aai=ql/V for' the small angles under con- 



sideration. The change in moment due to the change 
in the angle of attack at the tail is 



dC Xt 



AM^lpWSn^Aat 



where — =— slope of normal-force curve of toil surfaces. 
da, 

S„ area of horizontal tail surfaces. 
i)„ tail efficiency, 
or, in coefficient form, 

AC bt =- ^ 

Introducing an empirical factor K to allow for wing 
damping 

dC n K l S,dO Xt I 
-djc -& Vr da7 V 

Therefore (7c) may be rewritten as 

Dividing (8) through by k Y * and multiplying by t to 
make the expression nondimensional, 



1 P dO Zt S t 



'2k? K 



1 roc dC m 



$6 



2 P £fc r » da ^-^W 



771 

Letting ^ = M (reference 3) 

1 P ^dC z ,S t 1 lc dO m a cP0 

For convenience let 



0) 



(10) 



and 



1 P dO Zt S t 
Vt 2k? ~da~i S" =Wl * 

1 lc dC m _ m 

2 £7 ~d~a~~ m " 

cPO 

Tn t Tg+pm a Aa=i 2 -^ 

de , 



(11) 



Since 6={a-\-y) and g=^(see fig. 1), 2 and 6 may 

be replaced by a and y reducing the number of vari- 
ables in equations (7a), (7b), and (11) to three. These 
equations are rewritten as, 

-CjAt^CdAV- 



dC D . , dV' 



C»Ay-2C^V>- d -^Aa=-2r*g 

m < T (ft + W +maAa= K^ 2+ d?) 



(12) 



The foregoing equations of equiUbrium must each 
be satisfied at each instant (neglecting approximations 
assumed) of flight. It is obvious that each of the 
variables Ay, Aa, and AV will affect each of the others, 
and it is reasonable to assume that, if any one of them 
follows a regular scheme of variation with time, then 
the others will vary according to the same scheme. 
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Assume therefore that the variables change accord- 
ing to the exponential expressions: 

Ay=Ae*' 
Aa=Be™ 
AV'=Ce>-' 

(In mathematical terms is an integrating factor.) 
where 

A, B, and O are constants depending on the mag- 
nitude of the initial disturbance. 

e, the base of natural logarithms=2.71828. 

X, an arbitrary constant. 

t, time in seconds. 
It appears that if X is positive A7, Act, and AV will 
increase with time and therefore the motion will be 
unstable. If X is negative the departures from the 
condition of steady flight will decrease and the motion 
will be stable. 
Now if 

Ay=Ae u 

etc. 
dC n 



Substituting in (12) 
-OvA^-lO^ 1 



■Be xt =2T\Ce" 



aa 
fin 

0 D Ae^-20 Il Ce"-^Be^=--2rA\e^ b 
m,T(BXe x, + AXe*9+ Mm a Be x, = XAe^c 
Each of the equations may be divided by e x ' 

-G i A-2C d C~£B=2t\C a 

GoA-2C L G-^B=-2r\A b 

m,,THA+B)+iiin a B=(T\y(A+B) c 
From (14b) 

2 T \A+C D A-^B 



(13) 



(14) 



2(7, 



Also, since t and X do not appear except in the product 
tX, let X'=tX 

Substituting in (14a) 

-G^-%(2yA + G^^B)- d £ B 

Multiplying by C L 

- CJA^GdK'A- cja+gJ^b- g^b 

= 2 (X') 2 A+ X' CdA- X'^5 



or, since Gl+G d 2 —O b 2 

A[+C B t +ZC D \'+2{\')-} 

Therefore, 

a dCi, rj dOp , ^, dG L 

B C B 2 +3C D \'+2(\ , ) i 
Also from (14c) 

A[m,\'— (X / ) 2 ]=-B[(X , ) 2 -m e X'- / tmJ 
Therefore, 

A_ /xm tt +m c X / -(V) 2 
5 -m J X'+(X / ) i 

It therefore appears that equilibrium at all times is 
possible when, 

n dOh n dC^. /dO^ 

° D ~d~^ ^rfa tA ~d^_nma+m,y- (X') 2 (15) 



G B 2 +ZC D \'+2(\'y -m 5 X'+(X') s 



or 



_ mt y+(xO*](c7.g-^+V§) 



=[C^+3C! D X / +2(X0 s ][/iro«+ro e X'-(X / )T 

Expansion of this relationship and collection of terms 
give the biquadratic expression 

(\>Y+ (x') 3 [-m,+J(3C D +^)] 
+ (X') s [-m J |(3^+^)- M m a * 



— /xm a ~2-=0 



(16) 



For simplicity the biquadratic may be expressed as: 

_A{\y+B(\'y+c(\'y+D\'+E=o (17) 



where 
A=l 

a=[-m4(3C- 0+ ^)-,m a+ A(^-^ + ^)] 

Unfortunately there is no simple, direct method of 
solving biquadratics. It is possible, however, to factor 
a biquadratic into two quadratics, each of which is 
susceptible to direct solution. 
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The biquadratic 

(X')*+5(X') 3 +<7(AT+£*'+#=0 
may be replaced by the expression 

[(X') a +aiX'+6 1 ][(X') J +a s X'+6 i ]=0 
from which 



or 



It appears that 



*'--f±V(f)'- s - 

B=a 1 -\-a i 
P=a 1 6 a +a s 6 1 



(18) 



The general case may be worked out as follows: 
a^—B — di 

C=a l (JB-a l )+b 1 +b i 
B=a l b 1 - s r(B-a l )b l 
E 



b 2 = 



Therefore, 



^7=o ] (B-a 1 )+6 1 +^ 
E 



D=af+(B-a i )b 
Dropping the subscripts: 



(19) 



and 



VB-bD 
b 2 -E 



(20) 



These relations can be solved by plotting the curves 
of b against a. There are two intersections of these 
curves, in general, corresponding respectively to ai, 
bi, and to a 3) 6 3 . 
Also from (19) 



E=C-b-a(B-a) 

D=aC-ab-a?(B-a)+b(B-a) 

D-aO+a?B-<P 
b ~ ~B=2~a 

Substituting (21c) in (19a) 

'D-aC+cfB-a? 



(21) 



C=a(B-a)+(j 



B-2a 



) 



or 



+ E (D-aC+l'B-a?) (22) 

a 6 -a 6 (35) +a i (3B 2 +2G) 
-a 3 (B 3 +4BO)+a 3 (2B 2 C+BD+C 1 -4:E) 

-a(B 2 D+BO i -iEB)+BGD-D 3 -B 1 E=0 (23) 



By use of the foregoing relationships the coefficients 
of the quadratics may be determined with as high a 
degree of accuracy as desired by graphical means or 
by trial substitutions. 

At first glance it appears that a 6th-power equation 
such as (23) would be harder to solve than a 4th-power 
equation such as (17). In equation (17), however, the 
complex roots must be obtained; whereas in (23) it is 
necessary to solve only for the real values. Equation 
(23) is useful from a practical standpoint chiefly in 
obtaining accurate values of a* by making trial sub- 
stitutions from approximate values obtained from the 
expression for <h given on page 6. Because of the 
very small value of ch it is generally not necessary to 
include the terms in (23) that contain powers of a? 
higher than the third. 

Significance of X'. — As appears in equations (18), 
there are possible either 4 real values of X', 2 real 
and 1 pair of complex values, or 2 pairs of complex 
values. The values of B, G, and D in the normal 
flying range of conventional airplanes are always 
positive because of the signs and magnitudes of their 
constituent factors. It is obvious that no positive real 
value of X' can satisfy the biquadratic unless E is 
negative but that if E is negative there is such a 
solution for X'. A positive real value of X' signifies 
an aperiodic divergence. If — fim a is positive corre- 
sponding to static stability, then E is positive and the 
biquadratic expression indicates no possibility of an 
aperiodic divergence. 

The values of B, G, D, and E are, in general, such 
that the solution for X' gives two pairs of complex 
values. It can be shown by mathematical reasoning 
not essential to this treatment that an expression of 
Re*'"— where X'=f ±i^' and where K, f, and ^' 
are constants and i=V— 1 — can be replaced by an 
equivalent expression, K'ef' cos 5) where K' 

and 5 are new constants. Therefore Ay=A e vt ' may 
be replaced by Ay=A'e r ''' cos (ft'— S) and similarly 
for the other variables. It appears that the motion 
indicated by a complex root is therefore made up of 
sinusoidal variations of the angle of attack, angle of 
attitude, angle of the flight path, and velocity along 
the flight path and that the amplitudes of the oscilla- 
tions increase or decrease with time depending upon 
whether f is positive or negative. From equations 
(18) it appears that 
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where ft', f 2 ', fa', and fa' express the stability charac- 
teristics of the airplane. 

From the expression Ay=A'e t ' t ' cos {yf/'f — 5) the 

period of the oscillation is P'=^-, and the time for 

the oscillation to damp to one-half amplitude from any 

0 693 

instant chosen as the time of starting is T'= '—, — 

Since X' was set equal to Xr during the derivation, P' 

and T' are in units of time which are equal to t=^ 

seconds or the time in seconds necessary for the air- 
plane to travel the length I at the velocity V multi- 
plied by the relative-density factor p.. Therefore the 
final expression for the time in seconds for the oscilla- 
tion to decrease to one-half amplitude is 



„ —0.693 -0.693 
1= — z> — t= 



m 



f pSV 



-0.313 /wr 

=— p— Vg-ft 



(24) 



under standard conditions; and the period in seconds is 



2.83 rwz 



(25) 



under standard conditions. 
Substituting for a in equation (23): 



(2f) 8 + (2r') s (35) + (2f ')W+2C) 

+ (2f') 3 (B 3 +4BC) + QD'G&G+BD+C-lE) 

+ (2f ') (WD+BC-iEB) +BCD-D I -B i E=0 



From equation (21) 



,_ z?+2r<?+3G-y.g+4(r) 3 



(26) 



(27) 



Reference to equation (26) reveals that from the 
values of B, G, D, and 2? in the normal-flying ranges the 
coefficients of the terms in 2f will all be positive. 
Therefore a positive real value of f can exist only if 



{BCD— If 1 — WE) is negative. This fact was first 
pointed out by Routh (reference 4) by a somewhat 
different derivation and the factor is known as "Routh's 
Discriminant. ' ' 

There has been developed in reference 1 an approxi- 
mation for the case when D and E are small with 
respect to B and G by assuming that 



Then 



or 



and 



a l =B 

D=Bb 1 +a 1 C 
E=Cb 2 



O 



_D BE 



It follows that, 



B 

fl - 2 



fs_ 2\G C*J 



(28) 

(29) 
(30) 



, , IE UD BE\* 

These approximate expressions are generally used 
and are satisfactory for most cases encountered. 

Since B and C are large in the normal-flight range, 
values of f/ and fa' define a very heavily damped 
oscillation which is, in general, of short period. This 
oscillation is unimportant, except possibly in special 
cases above the stall; throughout the remainder of the 
report f and tf/ will be used without the subscripts 
to refer to the slightly damped phugoid oscillation 
previously defined by f 2 ' and fa'. 

Derivation of expression describing the sinusoidal 
motion. — The expression Ay=Ae > ' t =Ae K ' t ' can be 
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replaced by the equivalent expression Ay=A'e t ' 1 ' 
cos (ip't'—8) where A' and 5 depend upon the instant 
from whence time is taken as zero. If time is zero 
when Ay is at a point of maximum amplitude A'=Ay Q 
and 5=0 so that 

Ay=Ay 0 e t ' 1 ' cos ft (32) 

It follows that 

Aa=Ao^ n ' COS (^'i'-5i) (33) 

Ay =AV 0 'et'" cos Wt'-h) (34) 
A0=A0 o e r '" cos Wt'-h) (35) 



where Aao, A TV, and A8 0 are the magnitudes that the 
variables would have were they at a maximum at zero 
time and Si, 5j>, and 53 are phase angles of the variables 
with respect to the flight-path angle. 

These values of Ay, Aa, and AV may be substi- 
tuted in equations (12) and solutions made for f' 
and f. It will be found that the complete solution 
for £' and ^' given in the preceding section can be 

checkedinthismannerandthevaluesof^^j tanS i> 

and tan 5j are found at intermediate steps to be: 



A70 



Ol- 



dC D 
da 



Aoq 



tan 5i= 



^+g+4T) 



(3C7 D +4f) C0S 5l 



\° D da 0l 



— +r^)-^^{m^+ f m a -i(r'y-{m} 



r 5^' rS da J' 



dO L 



dC D 



Ay 0 



,dCp\ 
da J 



sin Sa- 



tan Sa= 



dOr. 
da 
dC D 

da 

N 



cos 5^ 



where 



t )(c D+ ^1 



(36) 



(37) 



(38) 



(39) 



N = -2 M m a -4r'm t -2(^) 1 +6(f') s +(2r'- 



The detailed derivation is not repeated because it 
is quite long and tedious and introduces no new con- 
cepts. 

Application of mathematical formulas. — The rela- 
tionships derived in the preceding section of this 
report make possible various treatments of the prob- 
lem of deterrnining stability. The period and time 
to damp to one-half amplitude may be determined" 

directly having given G L , G D > ^7 7[a'~ m <" 

— pM a from the equations (17), (24), (25), and either 
(26) and (27) or (30) and (31). 

An understanding of the underlying principles gov- 
erning stability can be had only from a consideration 
of the variations of the angle of attack, angle of atti- 
tude, angle of flight path, and velocity along the fright 



path relative to each other. Figures 2 to 14 are in- 
cluded to show the nature of these variations. Figures 
2 to 9 present the effects of changing —m Q and —fim a 
independently for two sets of typical values of 
dO dO 

G L , C D , -j^> and Figures 10 to 14 illustrate 

qualitatively the variations in phase relationships re- 
sulting from changes in -fim„ and serve as a basis 
for the discussion of their effects upon the stability 
of the airplane. The relationships (36) to (39) can 
be used to determine the characteristics of the motion 
for any particular design. 

The solution for the stability characteristics gives 
no direct indication of the effect of variation of the 
individualparameters upon the stability characteristics. 
By means of the relationship that follows directly 
from equation (22) 



m ff (c+2f')+Mma-[d+2f'c+4(f') z ] 

f-m^+2rc+4(i-Q 2 ]-^m g ( e +2r)+2r / ^+2r / c+4(r , ) 2 ll 



+ 



-m 4 +(c+4f') 
-m,+ (c+4n 



-m I [d+2f'c+40-') 2 ]-MTO a (e-l-2r)+2r'[d+2rc+4a-' 



(40) 
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where 

e=1.5<7 D 

charts may be plotted showing the variations of — m t 
with — fim a necessary to secure given values of £' 
and A number of these charts have been prepared 
covering the range of conditions likely to be encoun- 
tered in normal flight and are included in this report. 
(See figs. 15 to 54.) 

DISCUSSION 

The mathematical relationships evolved in the pre- 
ceding paragraphs permit calculation of the probable 
stability characteristics of a proposed design, but they 
offer little information as to the relative importance 
of various factors or as to the reasons for the effects 
produced by changes in those factors. In the following 
paragraphs the oscillatory motion is first considered 
in detail so that the sequence of flight conditions that 
must exist if instability is to arise may be pointed out- 
Next is given a general discussion of charts (figs. 15 
to 54) that show the effects on the stability character- 

dO 

istics of the six fundamental parameters C L , 0 D , 

-~, —p.m a , and — m ff . Finally, the effects of various 

physical characteristics of the airplane on its stability 
are considered in the light of the earlier discussion. 

The oscillatory motion. — When the terms in AV 
were ehminated by simultaneous solution of the equa- 
tions of equilibrium tangent and normal to the flight 
path, respectively, the following equality was found to 
exist: 



dC L 



q dCp . y, dOi, 
L da da _ itmg-\-m ll \' — (\') a 



da 



The terms to the left of the equality sign arise out of 
the necessity for equOibrium of forces; the terms to the 
right arise out of the necessity for equihbrium of 
moments. If a case be investigated with — nm a = 0 
dC D 



n dO L 
L/D-5 — 
aa 



■Or.- 



do 



-V^=-C r B s -3C7 D X / -2(X') a 

Q/CC 



or 



OT .+l(3a,+^)v + i(c^-ct§+cy)=o 

dC 

Since both C D and are positive below the stall, it 

follows that if — vcm a is zero the linear motion of the 
airplane is a heavily damped oscillation with the period 

and the time to damp to one-half amplitude 

rpr 0.693 

all in nondimensional units, or the motion is an aperi- 
odic convergence with the time to converge to one-half 
amplitude 



0.693 



depending on the relative magnitudes of the quantities. 

Since instability of the linear motion cannot arise 
when —iim a is zero, it seems probable that under- 
standing as to the underlying causes of instability may 
be gained by considering the effect of —fnn a upon 
the sequence of flight conditions during oscillatory 
motion. In the following paragraphs the part played 
by —nvia is analyzed by physical reasoning. Only 
the case when —nm a is positive is considered because 
when —iitn a is negative an aperiodic divergence from 
the equilibrium condition occurs and the question of 



dynamic stability does not arise. It is first pointed 
oat that —\im a plays a primary part in determining 
the phase angles between the change of angle of attack 
and the changes in attitude and flight-path angles. 
It is next shown that stability of the angular motion 
depends upon the phase angle between the angle-of- 
attack change and the attitude change. It is finally 
shown that instability of the linear motion can arise 
only when the phase angle between the angle-of-attack 
change and the flight-path-angle change falls within a 
certain range. 
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Fiookks 2 and 3.— Variations of components of longitudinal motion with time. 
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Figures 10 to 12 illustrate qualitatively the effect 
of -fii» a upon. 8i and 5 4 , the phase angles between the 
angle of attack and the flight path and attitude angles, 
respectively. When — pm a is zero the attitude does 
not change with time and hence, if an oscillation is set 
up in the flight-path angle (7), it is obvious that the 
angle of attack must change as shown in figure 10 where 
5i is 180°. Although not shown in figure 10, physical 
reasoning leads to the conclusion that 5 4 is 90°. When 
—firria is very large the moment tending to cause the 
attitude to change with the flight path is also large 
and the angle-of-attack change is small. The condi- 
tions are approximated by figure 12. The phase 
angles 5i and S t are also small and approach zero as 
—pm a is increased. Figure 11 shows an intermediate 
case which more truly represents the usual condition 
than either of the others. Here S^ISO 0 and 5 4 <90°. 
The curves of figures 10 to 12 were drawn from physical 
reasoning but are fully confirmed by the computed 
curves of figures 2 to 4 and 7 and 8. 

Referring to the curves of figures 2 to 9 one sees 
that, if 2 points are chosen exactly 1 period apart when 
A7 is zero, the value of A V will be found to be greater 
at the second point than at the first if the motion 
is unstable, to be unchanged if the motion is neutrally 
stable, and to be less if the motion is stable. Similar 
observations can be made with respect to the val- 
ues of A7 when AV is zero and to the values of 

^g-when A0 is zero. In other words, when neutral 

stability exists, the momentum along the average 
flight path m(V+AV), the momentum perpendicular 
to the average flight path (mVsin A7), and the angular 
momentum are each the same at the end of a period 
as at the beginning. If instability exists, there is 
more momentum at the end of the period for each of 
the types of motion than at the beginning. It is 
therefore desirable to consider the sequence of flight 
conditions that results in the increase or decrease of 
momentum in the system. 

Consider first the angular motion. Figure 13 is an 
illustration of a typical set of variations of Aa, A7, 
and Ad when f ' is zero. At the beginning of the half 

period when Ad is a maximum, Ag ^or A is zero and 

the angular momentum BAg is therefore zero. At 
the end of the half period Ag is again zero. It there- 
fore follows that the angular momentum introduced 
into the system during the half cycle must be balanced 
by the momentum removed during the same interval. 

At each instant when there is a finite value of Aa, 
there is a pitching moment introduced into the system 
1 dO 

equal to Aa-^pV 2 Sc-^- Since a moment M applied 

during a time t introduces angular momentum into the 
system equal to Mt 



=J o ±pV*Sc^Aa<fit' cos (ft-Bjdt 

represents the angular momentum introduced into the 
system during one-half cycle (neglecting the effects of 
terms involving the products of two small quantities, 
as has been done throughout the discussion). 

Considering, for simplicity, the case of neutral 
stability 

1 dO ~1 

ipV'ScAoo^f 

=- ; —sin (#-S«) 

pV^cAoo^^sin S t 



is the angular momentum introduced into the system. 

The rotational velocity causes a change of angle of 
attack at the tail that introduces a moment tending to 
oppose the rotation. The magnitude of this change is 

, dd -jdd 
. lAq_ lA Tt L Tt 



since 



.d$_de 

*dt~dt 



do 



) 



The moment is — Kr),^pV 2 S I ~g)~Y' ' ^ ee P- 3 - 

The angular momentum removed from the system 
Q"=0) during the half cycle is 

In the case of neutral stability the momentum intro- 
duced must equal that removed and therefore 



pV 2 ScAcJ^^h 



/dC Xt S\ 



or 



/im a Aoo . 

—TT 3111 

4> A0 O 



-m-=0 



(41) 



From figures 2, 3, 4, 7, and 8 it is evident that as 

— um a increases 4r> and 5 4 all decrease. The first 
4/ A70 • 

term in equation (41) may therefore either increase 
or decrease with increase in —pm a ; hence increasing 
— pm a may decrease the angular momentum added 
during a cycle. What happens in an actual case 
depends on the magnitudes of the quantities, as will be 
brought out later in the discussion. 
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Fiotoe 10.— Variations of angle of attaok, angle of pitch, and angle of flight path 
when static stability Is »ero. Damping coefficient assumed lero. 
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jFiqtjkb 13.— Variations of angle of attack, angle of pitch, and angle of flight path 
when static stability Is great. Damping coefficient assumed lero. 
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Fioube 11.— Variations of angle of attack, angle of pitch, and angle of flight path Fiowre 13.— Variations of angle of attack, angle of flight path, rotational velocity, 
when statlo stability is small. Damping coefficient assumed lero. and angle of pitch for an Illustrative case. Damping coefficient assumed zero. 
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Figure 14. — Variations of angle of attack, rate of change of angle of flight path, and 
angle of flight path for an illustrative case. Damping coefficient assumed lero. 
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The introduction of linear momentum along and 
perpendicular to the flight path might be discussed by 
considering equations (12a) and (12b) separately but 
the analysis may be simplified by combining the two 

dV 

through elimination of terms involving AV and 
The relationship established is: 

_^A 7 -3^-2^= 



dC, 



dC D 



V doc da , 



> 



dC L da. 
da dt 



Considering the case of zero damping during the in- 
terval from ypt=0 to yf/t=v it is obvious from figure 14 
that the integral sum of terms multiplied by A-y or by 

-~ is zero because these quantities are symmetrical 
dt 

with respect to the zero line. The equation of equality 
of momentum introduced to that removed becomes: 



* i ° T T dO iL dajt 
^i=o da dt 



which reduces to 



-ZO D 



da /Aao . . . dC L Aoq n ^ x n 

--r-^ sin 1 — cos o^U 

Ay 0 da A70 

(42) 



(' 



The growth of linear momentum therefore depends 

upon dij and 4>' in addition to the parameters 
Ayo 

O d , C L , and Since ^ is normally positive 

the following facts appear: G D is a factor serving 

always to remove momentum from the system; 

dO dO \ 1 
Gd-j^— Gl-^-j is a factor adding energy to the 

dG dG 

system when ^~J £ ^ > ^o~^> 05 is normally the case 

dG 

at high angles- of attack; and adds momentum to 

the system when 5!<90° and removes momentum 
when 5!>90°. When —fim a is zero, 5!= 180° and the 

factor removing energy from the system is (^^-f - 

as was shown earlier in the discussion. As — tim a is 

AcKn 1 

increased, hrj^t and all decrease (figs. 2, 3, 4, 7, 

and 8). The factor (cJ^-G^)^^ sin* 

reaches its maximum at some point between 5!=180 o 

and 5!=90°. The factor ^ ^ cos 5! changes from 

da A7 0 

a stabilizing to an unstabilizing influence at 5! =90° 
and reaches a maximum as an unstabilizing element 
between 5!=90° and 5j=0. It would appear at first 



dG 

thought that the effect of as an unstabilizing 

factor would always be small because of the decrease of 

^with increase of cos 5i but such is not necessarily 
A?o 

the case because of the large ratio of to G D , particu- 
larly at low angles of attack. 

There is an interesting point in connection with the 

(dO dC \ 
Gd-^—C l -^\ that appears as the chief 

unstabilizing element at small values of —y.m a . If the 
law of induced drag be applied to determine G D and 

rg^> ^ appears that 



( n dC L r dC D \_ p dG L , r dO L _ r dU E 

n ? dG L 
JL 1a 



_ r dC L 



"5 



■Or; 



dO L 
da 



2C L ^ L 



dC, 

~~d~a 



/ rt Cx, S \ 



Therefore equation (42) may be rewritten as 

~ z0d+ iH? ^} 0d ~ 0d ^W 5l+cos 5l ] = '° (43) 

It appears that instability of the linear motion for 
small values of —nm a can occur only at angles of 
attack where C D ^>G Dp . The detrimental effect of 
0 Df is not as great as at first appears, however, bo- 
cause the damping factor 3 G D increases with G D{ also. 

In the foregoing analysis, the effect of f ' has been 
neglected and therefore the statements made cannot 
be considered rigorously true when the stability or in- 
stability is of appreciable magnitude. The analysis 
does, however, point out the more important influences 
and the nature of the interaction of events that brings 
them into play. Equation (43) cannot be applied near 

dC L 



dG, 



the stall where ceases to equal 



s 



In figures 2 to 14, inclusive, and in the discussion of 
the oscillatory motion an arbitrary value of Ay 0 = 1 bas 
been assumed. No attempt will be made to evaluate 
A-y 0 in terms of control movements or gust velocities 
as the actual magnitude is unimportant with respect 
to stability (assuming that the deviations are not so 
great as to destroy the validity of the basic assump- 
tions). It seems not out of place, however, to suggest 

that relative magnitudes of ^ ^-7- S h etc., 
" A70 A70 "Vo 
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may have important bearing upon the comfort, ease of 
handling, and design load factor of aircraft and that 
further work taking these factors into consideration 
might lead to valuable information. 

Charts of rotational damping factor against static 
stability factor. — In the section dealing with applica- 
tion of mathematical formulas a convenient graphical 
means of showing the variation of f ' and ^' with — m t 
and —iim a has been described. This method of pres- 
entation was given by Gates (reference 5) in essen- 
tially the same fashion although differing in detail. 
In reference 5 preference is given to charts with 
coordinates of tail volume and fore-and-aft location 
of the center of gravity. Preference has been given 
in this report to charts with coordinates of rotational 
damping factor and static-stability factor because they 
are more convenient for use with data from wind- 
tunnel tests, they do not require assumption of arbi- 
trary fixed values for several important factors, and 
they cover a much wider range for a given number of 
charts. The charts are entirely nondimensional and 

vary only with 0 L , 0 D , -j^> and Variation of 

empirical factors such as the K in — m 1} ij t , ^ (see 

p. 17), etc., affect only — m„ or —iim a , as the case may 
be, and their effects upon the stability are readily 
apparent. 

The series of charts presented (figs. 15 to 54) are 
intended to show somewhat more precisely and com- 
pletely than has been done in the preceding discussion 

the effects of the parameters G L , 0 D , 

upon dynamic stability and to provide a convenient 
graphical means by which the designer may estimate 
the probable stability characteristics of a proposed 
airplane and the effects of various changes without 
recourse to extensive calculations. The charts cover 
the range of values of the parameters that appear 
likely to be attained in the near future. The values 
represented are summarized in table I. 

Certain general characteristics of the charts are im- 
mediately apparentuponinspection. As — Atm a increases 
from zero, f at first becomes more positive correspond- 
ing to a decrease in dynamic stability but, at a fairly 
small value of — fim», changes its trend and b ecomea more 
negative. This tendency is general throughout all the 
charts and is apparent whether — m q is large or small. 
It will be remembered that such an effect appeared 
probable from the discussion of the introduction of 
angular momentum into the system during a cycle 
by — m„. As pointed out at that stage of the report, 

sin 5 <( and ^> all decrease with increase of — iim a ; 

the charts of variation of the stability characteristics 
with — m, and —fim a show definitely that the decrease 
of the product of these factors is much more than 
sufficient to nullify the increases in - M m a after a 
certain value of —nm a has been exceeded. Increase 



of hence decrease in period of the oscillations, with 
increase of — /tm a occurs at practically all values of 
— m 4 and — iarn a . 

Increases of — m t give, in general, more negative 
values of f ' and consequently more rapid dying out of 
the oscillations. At large values of —iim a increases 
in — m, give more positive values of f but this effect 
is not of practical importance. It will be noticed that 
at fairly large values of — pm a increasing -m t has but 
slight effect but that -m, becomes of increasing im- 
portance as — jum a is made smaller. Increasing — m„ 
decreases rather gradually, giving oscillations of 
longer period. It appears that if the criterion of sta- 
bility be taken as the number of oscillations necessary 
for the amplitude to decrease to one-half its original 
value then the value of — m, is of particular importance. 

The charts of figures 23, 27, 28, and 29 show that 
increasing C L without changing other factors increases 
quite markedly the tendency to instabifity. This ob- 
servation agrees well with the effect to be expected 
from increasing Cl that might have been predicted from 

the part played by C L in the factor (cJ^-O^^ 

in equation (42). 

The effect of increasing C D without changing other 
factors appears in figures 29, 30, and 31. It will be 
seen that increasing 0 D reduces the range of instability 
and brings the curves of f'=& nearer together with a 
very large net increase in damping. The effect upon 
is negligible. 

The slope of the lift curve is of importance although 
not giving such extreme effects as changes in G L and G D 
(figs. 34, 35, and 36). Here again the effect is in good 

dC 

agreement with equation (42). Increasing tends 

to extend the instabifity region to greater values of 
—fim a or, in other words, if ~ixm a is large, 5 t is less than 
dO 

90° and becomes a definite unstabilizing factor. 

is a stabilizing factor, 



At small values of — m„ 



Increasing also 



as is apparent from the curves. 

tends to decrease 

Increasing the slope of the drag curve is distinctly 
unfavorable to stability as shown in figures 42, 44, and 
45. The period of the oscillation is but slightly 
affected. 

Reference to charts 15, 21, and 33 shows that the 
general effect of increasing angle of attack is to in- 
crease the region of instability but to decrease the 
spread between curves of £'=k. It appears that 
instability of the oscillatory motion is very unlikely 
for small angles of attack and becomes increasingly 
likely with increase of that angle. On the other hand, 
if — pm a and — m„ remain unchanged throughout the 
change of angle of attack, the value of f' may become 
more negative, which seems to be directly contra- 
dictory to the preceding statement. The explanation 
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lies in the fact that at small angles of attack changes 
of — nrn a and — m t have little effect upon the damping 
coefficient; whereas at large angles of attack a small 
change in — p.m a or — m q may change the damping 
coefficient from a definite negative value to a definite 
positive value. The period of the oscillation decreases 
markedly with the increase of angle of attack but the 
decrease is not so great as would appear to be the case 
from consideration of the increase in if/. The period 

in seconds is proportional to and the increase in 

0 L is sufficiently great partly to counteract the increase 
in y. 

The effects of various physical characteristics of 
the airplane on the stability. — In the preceding dis- 
cussion no consideration has been given to the various 
dimensional characteristics, aerodynamic interferences, 
etc., that determine the values of the fundamental 
parameters. A large number of factors affect the 
stability but in many cases the effects are of a minor 
nature. Only the more important ones will be dis- 
cussed in the following paragraphs. 

The wing loading appears in but two places in the 
analysis of stability in this report. The nondimen- 
sional relative-density factor 




I 



(standard conditions at sea level) 



is directly proportional to wing loading and the factor 



v 7? ^ L 

= o~oo — is proportional to the square root of 



wing loading. Since /x appears only as the coefficient 
of — m a it is apparent that the variation of the non- 
dimensional values of f 7 and \ff' with p are the same as 
for variations of —pm a . It therefore appears that, 
if the longitudinal motion of an airplane with respect 
to the air is to be unaffected, an increase in the wing 
loading of a given design must be accompanied by a 

dO 

proportional decrease in provided that I, c, and 

k Y are unchanged. Aside from the effect upon — fim a 
the wing loading affects the stability characteristics 
through the time factor t, which appears in the reduc- 
tion of the nondimensional units f and ^' to the time 
to damp to one-half amplitude and to the period of 
the oscillation in seconds 



-0.693 



and 



w 



If — p.m a is preserved unchanged as is increased, 

both the period and the time to damp are increased 
but the number of oscillations for a given degree of 
damping is unchanged. The effect of an increase of 
wing loading without change of other factors is there- 
fore to increase both ^' and r with the net effect, in 
general, of increasing the period; to make f more 
negative if — /im a is large, with a net effect upon the 
time to damp either of an increase or of a decrease; 
and to make f ' less negative if —y.m a is small, resulting 
in a comparatively large net increase in the time to 
damp to one-half amplitude. 

The aspect ratio of the wing combination affects the 
stability indirectly through its effect upon the para- 
ge dC 

meters and G D at a given value of O h . 

dO 

Increasing the aspect ratio increases and decreases 
dO 

both and 0 D . Comparison of the curves of figures 

36 and 39 corresponding to aspect ratios of 5 and 8, 
respectively, shows that (neglecting the effect due to 
dO 

change of -j^> which is comparatively small and will 

not materially affect the conclusions to be drawn) the 
net effect of increasing the aspect ratio is to increase 

the range of instability. The decrease in an 

effect favorable to stability, is less important than the 
decrease in 0 D , which quite markedly tends to decrease 
the stability. 

The effect of parasite drag is brought out quite 
clearly in figures 33, 36, and 40. (In general, increas- 
ing the parasite drag also increases the slope of the 
drag curve and this fact has been taken into account 
in the figures by using the method of reference 6.) 
It is apparent that the drag is an important item in 
determining the stability characteristics of a design 
and that the "cleaner" the airplane the greater the 
tendency to unstable oscillations. Increasing the 
parasite drag results in a decrease in the unstabilizing 
factor (Pd— 0 Dp ) and an increase in the stabilizing 
factor 3<7d, which appears in equation (42) in the 
discussion of the growth of linear momentum during a 
cycle. It appears therefore that air brakes of various 
kinds may be expected to have beneficial effects upon 
the longitudinal stability provided that they do not 
introduce undesirable moments or interferences. 
Highly efficient designs that achieve high values of 
Ginnix "without the use of devices introducing parasite 
drag may be expected to be deficient in longitudinal 
damping. 

The moment of inertia of the airplane about the 
lateral axis appears as the term k r 2 in the denominators 
of the expressions for -m, and —m a . Increasing lc Y 2 
decreases both -m t and -itm a proportionally with a 
corresponding increase in the time to damp. The 
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effect may be readily visualized by imagining any 
point -ism a , -m q upon the chart moved from or 
toward the origin along a radial line an amount such 
that the ratio of the new distance to the origin to the 
original distance is the ratio of the original value of 
k Y 2 to the new value. It is therefore to be expected 
that distribution of the mass of an airplane along its 
longitudinal axis may lead to unstable oscillations 
although other factors such as length and size of tail, 
location of the center of gravity, etc., are such that 
ordinarily the damping would be satisfactory. The 
effect upon the period of changing k T * is not great 
although increasing k Y * will, in general, lead to oscil- 
lations of shorter period. 

The location of the center of gravity affects the 
longitudinal stability through its effect upon the 

term-^ appearing in the expression for — m a . The 

pitching moment of an airplane about the center of 
gravity may be expressed as 

= - 0 Zu (0- O a ) + O m0w + G Xa \ + O np 

, dO Zt is t 
+oc '~da~ t 'S^ v ' 

Differentiating with respect to angle of attack 

dO u _ -dO Za dO XaZ dC„ p 

ia~ — d~*~ ( C '- G J +-^r ~ c ^~dir 

, da,dO St is t 



da da, c 
or, for practical purposes, 



dO m _dO L(n dC Xa2 dC mp 

, da t dO Zt lS t 

where 

0„ ratio of distance of c. g. from the leading edge 

of the wing to the wing chord. 
O a , ratio of the distance of the aerodynamic center of 

the wing from the leading edge of the wing to 

the wing chord. 
C Za , normal-force coefficient of wing. 
0 Xa , longitudinal-force coefficient of wing. 

2, distance of mean wing chord below c. g. 
O mj) , parasite pitching-moment coefficient of fuselage 

and landing gear. 

■p*, ratio of change of angle of attack at the tail to 
a change of angle of attack at the wing, an em- 
pirical factor depending upon downwash (refer- 
ence 7). 
dC 

Since is positive, increasing 0„ by moving the 

dO 

center of gravity toward the rear tends to make-^^ 
less negative and to decrease —fim a . The effect of 



longitudinal center-of-gravity location upon the sta- 
bility characteristics thus follows from the preceding 
discussion of the effect of changes in — y.m a . The 
location of the aerodynamic center of the wing, de- 
fined by O a , affects the stability in the same manner 
as the location of the center of gravity but with the 
opposite sign. The value of O a given by wing theory- 
is 0.25 (reference 8). Wind-tunnel data indicate the 
actual values to range from 0.23 to 0.25 for conven- 
tional airfoils (reference 9). 

The vertical location of the center of gravity with 
respect to the mean aerodynamic wing chord is not 
unimportant in many cases and should be considered 
particularly in high-wing or low-wing monoplanes. 

Since is positive, positive values of (as in the 

dO 

case of a low-wing monoplane) make less negative 

and may cause an unstable divergence (static insta- 
bility) in cases that would be considered stable on the 
basis of calculations involving only O a . The value of 
dO x „. 

—j — is derived as follows: 
da 



0 Xw = — C Da cos a+0 L sin a 
da V 



/j> «A i ( n i dO L \ . 

-jcos -rf^jma 



dO D 



where 0 Bw and > *° are for the wing alone. For a case 

An is 
with 61=1.4, ^=4, (7^=0.01,^=6, «=15° 

67^=0.01+^=0.11 
dO D „ 2X1.4X4 



da 



"6?" 



-=0.60 



dO x . 



^=(0.8X0.966)4- (4.11X0.259) 
0.774-1.06=1.83 
Assuming the case of a low-wing monoplane with -= 

0 

0.25 and C„=0.30, it will be found that the stability 
will be the same as though 

^=0 and <7„= 0.30 +(i|^X 0.25^ =0.41. 

The example represents perhaps an extreme case but 
illustrates the desirabffity of investigating the effect of 
vertical center-of-gravity location. 

The value of the slope of the curve of parasite 
pitching-moment coefficient against angle of attack 
may be determined with reasonable accuracy only by 
careful testing of a scale model or the complete air- 
plane. For a number of designs of military airplanes 
for which sample computations were made and com- 

dO dO 
pared with measured values of — values of ," p 

da da 

were found to range from 0 to 0.4, with an average 

value of 0.2, at angles of attack corresponding to 

dO 

cruising and high speeds. The value of . mp became 

da 
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less positive as the angle of attack was increased in all 
cases where it was of appreciable magnitude but in 
no case became negative. It therefore appears that 
the parasite pitching moment will be, in general, such 
as to decrease — ftm a with corresponding effects upon 
the stability characteristics. 

The tail length, size, and aspect ratio are very 
important factors in the stability characteristics. Tail 
length appears both in — m e and in — ftm a . It will be 
remembered that 

P S, dO Zt 
da. 



^ky 2 S 

_ m Ic VdC Lfn n . ,dC x 2 

da, dCz, Sj 
^ITa dat S c 



da 



Increasing I has a very beneficial effect upon the factor 
— m e , which is proportional to the square of I. In 
an actual case increasing I will also increase k T , and 
perhaps in some cases rj,, so that the net increase in 
-m, will only approximate that indicated by consider- 
ing only the square of the tail length. Increasing I 
makes — p.m a larger because of the increase in moment 
arm of the tail. There are also secondary effects because 



I affects tt^j ij„ and k T . 



The net effect of small in- 



creases in tail length can only be predicted from con- 
sideration of the particular design. There is a small 
range of conditions, corresponding to small values of 
—fitjia, in which increasing I increases the time to damp 
because the effect upon —fim a is of more importance 
than the increase in — m q . Large increases in I will 
always increase the period and decrease the time to 
damp. 

It has been mentioned that I affects rj t and 

The effect upon rj, is small and may be neglected. 

da 

Increasing tail length increases -jr- The angle of 

attack at the tail may be expressed as 
a,=a-\-i, — e 

where i, is the incidence of the tail referred to the wing 
chord. 

The downwash may be expressed as 

«=^(x+l)-°- 38 (2/+l)-°- 23 C r i (reference 10) 

where R is the equivalent monoplane aspect ratio of 
the wing combination and x and y are the distances of 
the tail plane behind and above or below the trailing 
edge of the wing in chord lengths. The value e there- 
fore decreases exponentially with tail length and hence 

da, ^ de 

da d~a 

aa 

becomes greater with increase of I. 
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Fioubes IS to 20.— Dynamic stability charts. 
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Figtjbes 21 to 28.— Dynamic stability charts . 
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Figures 27 to 82.— Dynamic stability chart*. 
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Fiqubes 45 to 50.— Dynamic stability charts. 



--2.50 




O 10 20 30 

Static, stability factor, -lim^ 



<7l=2.40 
dCt. 



namz 49. 



Cd^O&S 
dCo 



=2.00 



6 

«s 

t 

0 

.0 J 



















-A 












-S' 






1.0 










-r- 




























F 












i 
i 

—i 








I.S 




1 




/ 






/ 

/ 

/ 

/ 






1 
l 


1 
1 




i 










1- 

/ 


1 
/ 














li 
1 

/ 1 


A 

' \ _ 
; \ 


.3 












i 
l 
i 


/ \ 

1 

1 


/ 
/ 










1.7 


1 


J-.2 

\ / 














/' 


— V 

A 
/ \ 














iff 




































\\\ 











10 20 30 

Static stability factor, -fim a 

Ci-2.40 Co =1.00 

<e*=«0 d -£~2 60 



40 



da 



nanse so. 



ANALYSIS OF LONGITUDINAL STABILITY IN POWER-OFF FLIGHT 313 




nooBK si. nouns is. 




314 



REPORT NATIONAL ADVISORY 



COMMITTEE FOB AERONAUTICS 



The size of the tail appears in the ratio SJS occurring 
in the expressions for -m, or— m a . Increasing SJS 
increases both — m t and — ian a . In the case of -m, 
the increase is proportional to SJS but in the case 
of —{im a the increase is more than proportional to 
the increase in SJS because — /im a would normally be 
negative were SJS equal to zero. Increasing tail size 
may be expected to have a net effect somewhat less 
beneficial than increasing I unless the airplane already 
has a large value of — ixra a . 

The aspect ratio of the tail surfaces affects the 

stability through Increases in ha, ve the s ame 

effect as increases in SJS and need no further discus- 
sion. It will be noticed that if SJS is increased by 
dC z . 

increasing the chord, —r— 1 will be decreased so that the 
da t 

net effect will be less than would be expected on the 
basis of tail size alone. 

The tail efficiency factor i\ t affects the stability 
characteristics in the same manner as the ratio SJS. 
Tail efficiency depends upon the interference effects of 
the wings, landing gear, and fuselage upon the tail 
plane. Warner states (reference 11) that, in general, 
tail efficiency has been found to range from 0.6 to 0.8 
with extreme values as low as 0.54 and as high as 1.2. 
Computations made by the author for a number of 
military airplanes indicate that for modern designs the 
tail efficiency factor will be of the order of 0.75 to 0.8. 

The series of charts of stability characteristics against 
— m 4 and —/wia has been extended to include the 
range of high-lift devices known at the present time 
(figs. 41 to 54). The charts bear no evidence of in- 
creased difficulty with longitudinal stability by the 
use of these devices. Devices which increase the 
parasite drag as well as the maximum lift coefficient 
are more desirable than those which merely serve to 
maintain the flow to high angles of attack. It will 
be noticed, however, that in order to use the charts 
with variable-area lifting surfaces the true area of the 
extended wing must be used in order to come within 
the range of parameters given. When the true area 
is used, SJS and /t will have been decreased and the 
quantity c(C 0 —G a ) changed. The effective value of 
-m, is thereby decreased as compared with the wing 
with lifting surfaces contracted; whereas the value of 
—ixm a may be either decreased or increased, depending 
upon the relative values of the quantities. The fact 
that an airplane has good stability characteristics 
with conventional surfaces is therefore no assurance 
that trouble will not be experienced when high-lift 
devices, particularly those increasing the wing area, 
are added. In connection with calculations of stability 
with high-lift devices it should be noted that in cases 
of such devices as split flaps, variable-area wings, etc., 
which greatly increase the camber of the wings, par- 
ticular attention should be paid to reductions in tail 
efficiency ij, which may occur in some instances 



because of the change of effective location of the tail 
surfaces relative to the wing chord. (See also refer- 
ence 12.) 

The analysis and discussion have been confined to 
consideration of only the long-period oscillations or 
the divergence of an airplane as a rigid body. The 
effects of such physical characteristics as wing or 
tan-supporting structure elasticity or of free longitudi- 
nal controls are not, however, unimportant. No 
attempt will be made to analyze these effects mathe- 
matically but they will be briefly discussed in connec- 
tion with the basic relationships brought out in the 
foregoing portions of the report. 

The effects of wing elasticity have been treated mathe- 
matically in reference 13. .Although these effects are 
normally not important they may in some oases, 
particularly in view of the high speeds being obtained 
with modern aircraft, be sufficient to change a nor- 
mally positive value of -/im„ to a negative value with 
possibly serious tendencies to divergence. The obvi- 
ous means of guarding against such possibilities are 
to make the wings rigid in torsion and to use airfoil 
sections having positive or zero values of C^,. 

The effect of elasticity of the tail-supporting struc- 
ture is primarily of importance in connection with the 
phenomena of buffeting and tail flutter and needs no 
consideration here other than to call attention to the 
possibility of resonance between the natural frequency 
of oscillation of the tail structure and the short-period 
oscillations of the airplane, which are normally very 
heavily damped. The results of such resonance would 
probably be structural failure and would not receive 
practical classification as instability. 

The effect of freedom of the longitudinal controls 
has been treated mathematically in reference 14. The 
mathematical analysis introduces several new factors 
into an already complex problem and it seems likely 
that no designer would attempt to predict the stability 
with free controls. From a practical standpoint it is 
obvious that if the control surfaces do not deflect 
appreciably as the airplane follows the sinusoidal vari- 
ation of angle of attack, etc., the stability will be the 
same with free controls as with fixed controls. 

The influences that cause the control surfaces to 
deflect are the aerodynamic hinge moment caused by 
change of angle of attack of the tail and air speed; the 
inertia moment that tends to keep the surfaces from 
following the changes of attitude of the airplane; and 
the mass moment, in the case of statically unbalanced 
controls, that tends to deflect the elevator under the 
influence of gravity and accelerations of the tail plane 
is accelerated normal to the longitudinal axis. The 
aerodynamic hinge moment, which is normally of such 

sign as to decrease and hence both — m ( and — /tm„, 

may be kept small by aerodynamic balance of the 
elevator surfaces. The inertia moment is small and 
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may be neglected, as in reference 14. The mass 
moment may be kept small by statically balancing 
the elevator. 

It follows from the foregoing discussion that sta- 
bility with free controls can be insured by providing 
ample stability with fixed controls and then aero- 
dynamically and statically balancing the elevator sur- 
faces or making the control irreversible. It is generally 
not desirable to attain complete aerodynamic balance, 
but the difference between the stability with free con- 
trols and that with fixed controls can be made small 
without complete aerodynamic balance or without 
making the controls completely irreversible. 

CONCLUSIONS 

1. It is believed that for the average student of the 
subject this substitute treatment of the derivation of 
the relationships governing longitudinal stability will 
bo more understandable than the classical treatment. 

2. Increase in the amplitude of longitudinal oscilla- 
tions with time is possible only when the change in 
angle of attack lags the change in flight-path angle. 
The lag of the angle of attack with respect to the 
flight-path angle depends primarily upon the static- 
stability factor —pm a . 

3. The nondimensional parameter — pm a not only 
plays the most important part in determining the lag 
of the angle of attack with respect to the flight-path 
angle but also is the only source of unstabilizing 
angular momentum with power off. Increasing tte 
value of —ftm a , i. e., increasing the static stability, 
will increase the tendency to instability if the original 
value of —nm a is small. After a certain value of 
— IstUa is reached, further increases in the negative 
sense will reduce the tendency to instability. 

4. The effect of the vertical location of the wing 
relative to the center of gravity should not be neglected 
in calculating the static stability. 

5. AerodynamicaJly efficient airplanes are more 
likely to be dynamically unstable than less efficient 
types. 

6. Increasing the wing loading will increase the 
period of the oscillations but will not necessarily make 
the airplane either more or less stable. 

7. High-lift devices that depend partly upon in- 
creases in wing area may lead to instability if the size 
and disposition of the tail surfaces are determined on 
the basis of the contracted wing. 

8. There is no likelihood of dynamic instability at 
the angle of attack corresponding to maximum or 
cruising speeds. The damping of the oscillations may 
be more or less pronounced at low speeds than at high 
speeds. There is likelihood of dynamic instability at 
low speeds. 

9. Stability with free controls can be brought close 
to that with fixed controls by aerodynamically and 
statically balancing the elevators. 



IL ESTIMATION OF STABILITY 
METHOD 

The charts included in this report (figs. 15 to 54) 
provide an easy method of estimating the dynamic 
stability of airplanes. In order to use the charts it is 
necessary to know the values of the nondimensional 

m 4 , and —pm a for the 

airplane at the speed under consideration. Knowing 
the parameters, turn to the charts and select those for 
dC dO 

which G L , C D) -j—— and -r^ most nearly represent the 

act da 

case under consideration. Spot the particular point 
representing — m c and -/im„ upon the charts. Eough 
interpolations between the values of f ' and ty' so in- 
dicated will, in general, provide a sufficiently accurate 
estimation of the stability characteristics of the design. 
The time to damp to one-half amplitude, in seconds, 
is 



parameters 0 L , C D , ^^>-^~ } 



-0.313 



r 



and the period in seconds is 

V 

The number of oscillations for the amplitude to die to 
one-half the original value will be p= — ^-p- 

DETERMINATION OP NONDIMENSIONAL 
PARAMETERS 

The nondimensional parameters should be taken 
from available flight or wind-tunnel data if possible. 
If such data cannot be obtained they may be calculated 
by methods outlined in the following paragraphs. 

Lift coefficient C £ - — The lift coefficient C L follows 
directly from the wing loading and the speed. In gen- 
eral, it is desirable to choose G L as 0.20, 0.80, 1.40, 
1.90, or 2.40 to facilitate use of the charts. For air- 
planes with variable-area wings the wing area in use 
under the flight condition should be used as the basis 
of calculations. 

Drag coefficient G D . — The drag coefficient C D may 
be estimated from the relationship 



Cd — Oj 



Dp- 



C 1 

(see reference 6) 



r s e 



where 0 Dp is the parasite-drag coefficient, 
(kbY 

^-g-> the equivalent monoplane aspect ratio 

(reference 7 gives k in convenient form), 
e, an airplane efficiency factor ranging from 1 
for very clean designs to 0.7 for very 
inefficient designs (reference 6) . 



316 



REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



dO 

Slope of lift curve — The slope of the lift 
dO 

curve may be taken as 4 (all values of angles are 

in radian measure) at low and intermediate angles of 
attack and as 3 near the stall if the values of the other 
parameters are based on the true wing area and the 
airplane is a conventional type. It should be noted 
that in cases of extremely high or extremely low aspect 
dC L 



^— differs considerably from 4 the 



ratios for which 

wing area may be considered such as to make 

dC 

x, Cd, and —j— i — m c , n, and — m a are 



^=4 if a 



da 



each multiplied by the ratio of 4 to the actual value of 
dC L 



da 



for reference to the charts. In such a case the 



ratio must not be neglected in the conversion to time 
in seconds. 



dO 

Slope of the drag curve 



-The slope of the 



drag curve may be estimated from the relationship 



dC D 



da 



S 



,da 



The value of -j— here used may be taken as 4 through- 
out the flight range as does not fall off with 

near the stall. 

Rotational damping factor -m t . — This parameter 
will, in general, have to be estimated although it may 
be determined from tests of an oscillating model, or 
from tests on a whirling arm (reference 15). For 
purposes of estimation the equation 

v P S t dO Lt 
-m q =K Vt 2£7^^; 

where K is taken as 1.25 will generally be sufficiently 
accurate. The tail efficiency t\ t may be taken as 0.75 
or 0.80 for modern designs. The distance from the 
center of gravity to the elevator quarter-chord point 
I, the ratio of tail-plane area to wing area SJS, and 
the radius of gyration of the airplane about the span 
axis k r are dimensional characteristics of the airplane. 
The radius of gyration may be found from the rela- 
tionship 



where 



B is the moment of inertia about the span axis. 
g, the acceleration of gravity. 
W, weight of the airplane, 
or estimated from the relationship 



where 

C B is an empirical constant. 
k, over-all length of the airplane. 
h, over-all height of the airplane. 
From 11 airplanes for which C B was determined 
(reference 16) 
0.0325^ C B ^ 0.0394 with the average value 0.0362. 
In the absence of test data the slope of the lift 
curve for the tail may be taken as 

~3^~ 5 *2 ( re ^ erence 17) 
1+ b? 

Static stability factor —pm a . — The static stability 
factor is 

W Ic dO nc ^. 

V™"- S P X32.2Xi 2E? da 

or 



W 

6.5^c dG m ,^ 



(standard conditions) 



ry- da 

where c is the chord length upon which 0„ is based. 
Various methods have been proposed for estimating 



dC, 



''■ The author prefers the relationship 



I S t dC Lt da t 
Vt c S da i da 



dC m . t dC Lf „ _ dC Xo z dO m „ 
da ~da KV ' ° a) + da c~ ] ~~3a~~ 
where 

O a is the distance of the center of gravity of the 
airplane from the leading edge of the reference 
chord in chord lengths. 

C a , the distance of the aerodynamic center of the 
wings (reference. 8) from the leading edge of 
the reference chord in chord lengths. 

dO Xw 

^ a > the slope of the curve of longitudinal force co- 
efficient for the wings alone against angle of 
attack and is 

dO Xa / dO Dw \ / dO L \ . 

la-^—da-) 003 tt+ [°^+-da-) Sm a 

where and 0 Da are values for the wing 
alone. 

— » the distance of the reference chord below the 
c 

center of gravity in chord lengths. 



dO, 



-j^> the slope of parasite pitching moment due to 

fuselage, engine nacelles, landing gear, etc., 
against angle of attack. 

the rate of change of the angle of attack at the 

tail with change of angle of attack of the air- 
plane and is 

j^=l — ^ where « is the angle of downwash at the tail. 
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The value of G a varies from 0.23 to 0.25 referred 
to the mean aerodynamic chord of conventional air- 

dC 

foils (reference 9). The value of * " may be esti- 
mated if the designer so desires. The author has 
found the value to vary from 0 to 0.4 for a few military 
airplanes at low angles .of attack with the average 

dO„ 

value about 0.2. The value of v " was found, for 

da 

the cases considered, to fall off to approximately zero 
at high angles of attack. 

The rate of change of downwash angle at the tail 
with change of angle of attack may be estimated from 
the relationship 

&=wr> ix+1) ^ (y+1) ^ (reference 10) 

where 
dO L 



da 



is in radian measure 



x, the distance of the tail plane to the rear of the 
trailing edge of the wing in chord lengths. 

y, the distance of the tail plane above or below 
the trailing edge of the wing in chord lengths. 

■n , » . _ de 0.65 dO L 

For the average case y=0, 3=2.5, < k e == 7j£j)2 ~2^~ 



EXAMPLES 
AIRPLANE A 



O,=0.80 



(7^=0.080 



dO, 



da 

dC D 
da 



^=3.95 



=0.39 



(Taken from full-scale flight data.) 

— m 4 =2.6 — jim a =16.5 

(Estimated.) 

From figure 23 f'=— 0.044 f'=0.48 
From figure 25 f'=— 0.035 ^'=0.48 
From interpolation between the figures on the basis 
of C D , £'=— 0.039 and ^'=0.48. From interpolation 

dC 

between the figures on the basis of f'= — 0.038 

and ^=0.48. From full-scale unpublished flight tests 
of this airplane ?'=— 0.041 and ^'=0.72. 



AIRPLANE B 



(7 £ =0.80 

<7 D =0.076 



dOi 
da 

dO D 



^=3.60 
=0.42 



-m t =1.5 



— nm a =9 

(Taken from full-scale flight data.) 
From figure 23 f'= —0.034 ^'=0.49 
From figure 24 f'= —0.030 ^'=0.50 
From figure 25 f =— 0.025 ^'=0.48 
From interpolation between figures 23 and 24 on the 
basis of C D , f'=— 0.031 and ^=0.50. From inter- 
polation between figures 23 and 24 on the basis of 
fin 

f =— 0.033 and \fr'=0.49. From interpolation 

dC 

between figures 24 and 25 on the basis of 

f =-0.027 and ^'=0.49. 

From full-scale flight tests of this airplane (refer- 
ence 2) f'=— 0.038 and ^'=0.55. 



AIRPLANE C 



Cr=0.80 



C»= 0.094 



dC> 



da 
dO D 



2=4.00 



=0.51 



-m 5 =2.3 



"3a 
— jim a =4.8 
(Estimated data.) 

From figure 23, f'=— 0.034 and ^'=0.4. 

No interpolation is necessary. From extrapolation of 

full-scale unpublished flight data f'=— 0.030 and \fr'= 

0.72. 



AIRPLANE D 



(7^=0.80 



Cd=0.12 



^=4.00 
da 



dO D 
da 



=0.54 



— /xm a =3.40 



— m e =1.76 
(Estimated data.) 

From figure 26, f =— 0.03 to —0.04 and ^'=0.39. 
From full-scale unpublished data f'=— 0.035 and ^'= 
0.53. 

Langlet Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Taangley Field, Va., December IS, 19S4- 



APPENDIX 



RECAPITULATION OF SYMBOLS 



A,) 
B, 

C, ^coefficients of stability quartic. 

D, \ 

j'jcoefficients of stability quadratic. 

e, base of natural logarithms. 

e, airplane efficiency factor sucb that 0 D t = 



C B , coefficient of resultant force. 

C z , coefficient of force normal to the wing chord, 
positive downward. 

Ox, coefficient of force parallel to the wing chord, 
positive forward. - 

C a , coefficient of aerodynamic center (ratio of dis- 
tance of aerodynamic center from leading edge 
to chord length). 

0„ coefficient of center of gravity (ratio of distance 
of center of gravity from leading edge to chord 
length). 

G mp , coefficient of pitching moment other than that 
from wings and horizontal tail surfaces. 

OmQ, coefficient of pitching moment of wing at zero' 
lift. 

Cm e .„ coefficient of pitching moment with respect to 

center of gravity. 
q _ ky 1 coefficient of radius of gyration about 
B ~W+&)' lateral axis. 
li, over-all length of airplane. 
h, over-all height of airplane. 
I, distance from center of gravity to quarter-chord 
point of horizontal tail surfaces. 

K, empirical factor by which computed value of 
rotational damping factor of tail is multiplied 
' to give total rotational damping factor. 
k, empirical factor for conversion of biplane span 
to equivalent monoplane span. 
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m 



lc dO, 



"e.g. 



static-stability factor. 



— m. 



2F" 5 ~S ~d~a' ro * a ^ ona ^ damping factor. 
0, subscript denoting that the value of the symbol 
refers to zero time. 
P' , period of an oscillation, nondimensional units. 
P, period of an oscillation, seconds. 
T', time' for an oscillation to decrease to one-half 

amplitude, nondimensional units. 
T, time for an oscillation to decrease to one-half 

amplitude, seconds. 
t, subscript denoting that the symbol refers to the 
horizontal tail surfaces. 
Ay, ratio of AV to V. 

w, subscript denoting that the symbol refers to the 
wing alone. 

" x, distance of horizontal tail surfaces to the rear 
of the trailing edge of the wing, chord lengths. 
y, distance of horizontal tail surfaces above or 
below the wing chord (extended), chord 
lengths. 

z, distance of wing chord below center of gravity. 
5i, angle of lag of change in angle of attack with 

respect to change in angle of flight path. 
5 S , angle of lag of change in velocity along flight 

path with respect to change in angle of flight 

path. 

S 3 , angle of lag of change in angle of pitch with 

respect to change in angle of flight path. 
Si, angle of lag of change of angle of attack with 

respect to change in angle of pitch. 
77,, tail efficiency. 
0.693 



f'=- 



damping coefficient. 



2tt 



f=jp' period coefficient. 

fi=-^j relative density. 

T= ^y' time conversion factor. 
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